Self-Consistent Modification To The Electron Density Of States Due To 
Electron-Phonon Coupling In Metals 
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The "standard" theory of a normal metal consists of an effective electron band which interacts 
with phonons and impurities. The effects due to the electron-phonon interaction are often delineated 
within the Migdal approximation; the properties of many simple metals are reasonably well described 
with such a description. On the other hand, if the electron-phonon interaction is sufficiently strong, 
a polaron approach is more appropriate. The purpose of this paper is to examine to what degree the 
Migdal approximation is self-consistent, as the coupling strength increases. We find that changes 
in the electron density of states become significant for very large values of the coupling strength; 
however, there is no critical value, nor even a crossover regime where the Migdal approximation 
has become inconsistent. Moreover, the extent to which the electron band collapses is strongly 
dependent on the detailed characteristics of the phonon spectrum. 

PACS numbers: 71.38.-k, 71.10.-w 



INTRODUCTION 

The Migdal approximation for the electron self energy 
due to the electron-phonon interaction consists of ne- 
glecting vertex corrections. This procedure was first jus- 
tified by Migdal [ij based on an approximate treatment 
of the first-order vertex correction. He found that the 
correction to the bare vertex is of order O(A^), where 
A is the dimensionless electron-phonon coupling constant, 
u>D is the typical phonon frequency, and ep is the elec- 
tron Fermi energy. If we ignore the factor of A Q , then 
the ratio ujd/^f is generally very small in a metal. 

Subsequently, Engelsberg and Schrieffer [3| performed 
numerical calculations of the self energy and spectral 
function, based on the Migdal approximation The 
result is found in several reviews and texts 0, , and we 
quote here the main results. The electron self energy is 
given by a frequency dependent, momentum independent 
function, 
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where ip( x ) is the digamma function and the entire 
expression has been written for a frequency just above 
the real axis, uj + i5. In Eq. (|TJ n(v) is the bose function, 
and a 2 F(v) is the electron-phonon spectral function. A 
truly self-consistent approach would require, amongst 
other things, a self-consistent correction to the phonon 
spectrum, due to the interaction with electrons. Migdal 
estimated this correction, and found that the phonon fre- 
quencies are renormalized, and an instability is encoun- 
tered as the bare coupling strength increases. However, 
we are adopting a more phenomenological approach here. 
The common practice || is to take information concern- 
ing the phonons from experiment as input into the theory 



for the electrons. The justification for this comes from 
experiment, where well-defined phonons are observed in 
neutron scattering experiments 0], for example. Since 
these are used in the theory for the electron properties, it 
would be incorrect to compute renormalizations for the 
phonons. We follow this philosophy in everything that 
follows [13 ■ 

Nonetheless, Eq. (JTJ was obtained with a number 
of other simplifying assumptions and approximations. 
In particular the self energy of the electron is deter- 
mined by an infinite set of diagrams in which phonon 
lines do not cross; these can be summarized by the dia- 
gram in Fig. 1, where the full electron Green function, 
G(k, uj + iS) = l/(u> + iS— £k — S(k, u> + iS)), is required, 
and of course depends on the very self energy that we 
are trying to calculate. Yet Eq. shows no sign of 
self-consistency. The reason was noted already in Ref. 
Q and arises because the bandwidth is assumed to be 
essentially infinite compared to the typical phonon en- 
ergy. Then, the nested diagrams which arise from iter- 
ating the equation in Fig. 1 all contribute zero, and the 
same result is obtained by simply replacing the full elec- 
tron Green function in the figure with the non-interacting 
Green function. 

Engelsberg and Schrieffer Q relaxed the assumption of 
infinite bandwidth ^l], but used realistic values for the 
phonon frequency and Fermi energy for materials known 
at that time. They found only very small effects. More 
recently, Alexandrov et al. [l3| readdressed the question 
of the impact of a finite (i.e. not infinite) bandwidth 
on the electron properties, and adopted much more ex- 
treme values of the ratio of the typical phonon frequency 
to Fermi energy, cod/^f (referred to hereafter as the fre- 
quency ratio). They concluded that the Migdal approx- 
imation breaks down for coupling strengths that exceed 
unity. 

In this paper we wish to assess this conclusion, by ex- 
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FIG. 1: Diagram for the electron self energy. Note that the 
self-consistent electron Green function (heavy solid line) is 
used in this calculation. 



THE SELF-CONSISTENT MIGDAL 
APPROXIMATION 

The self-consistent Migdal approximation results in the 
following equation for the electron self energy S(cj + iS): 



E(oj + i8) = / dva 2 F(v) 



-00 N o (0) 



n(y)+f(-u') , n(u)+f{u') 



lu + ic 



where 



N(u) = / deN (e)A(e,uj). 



(2) 



(3) 



In these equations f(oj) and n{v) are the Fermi and Bose 
distribution functions, respectively, N {e) is the non- 
interacting (bare) electron density of states, and N(lu) is 
the self-consistently calculated electron density of states. 
The electron spectral function, A(e, lu) is given by 



amining the effect of a more realistic phonon spectral 
function. Alexandrov et al. |l2j used an Einstein spec- 
trum to simplify the calculation. This spectrum is, of 
course, singular, and it is perhaps not too surprising if 
singular behaviour in the electron properties results. We 
will first outline the problem as posed by Alexandrov et 
al. and demonstrate that singular behaviour exists for 
any coupling strength. We argue that this behaviour 
does not necessarily invalidate the calculation, as tac- 
itly assumed in Ref. 01 • Instead those authors properly 
focused on the more global behaviour of the electron den- 
sity of states (EDOS), as a function of increased coupling 
strength. We examine in much more detail this global be- 
haviour as a function of the parameters in the problem, 
so that a more quantitative assessment of the breakdown 
can be obtained. In particular we examine the depen- 
dency of the band collapse on the frequency ratio, the 
electron-phonon coupling strength, the presence of a sec- 
ondary band, the shape of the bare band (Lorentzian ver- 
sus square), and finally, the shape of the electron phonon 
spectral function, a 2 F{y). For convenience a 2 F{y) is 
modified from an Einstein spectrum to a Lorentzian. In 
this way a single parameter (the width of the spectrum) 
controls its shape. We find that there is no clear transi- 
tion or even crossover to a regime where the Migdal ap- 
proximation has become inconsistent. Nonetheless, this 
conclusion is not meant to imply that this calculation 
shows that the Migdal approximation is accurate in the 
intermediate or strong coupling regime. As will be sum- 
marized in the final section, other work suggests that 
this is not the case. Our calculation merely shows that 
within the Migdal framework, a signal of this potential 
breakdown does not occur, if a broad phonon spectrum 
is used. 



A(e,w) = lmG(e,uj + iS), (4) 

where the single electron Green function is given by 



G{e,uj + id) 



lu + iS — e — Yj(w + id) 



(5) 



Notice that we have tacitly assumed that the electron self 
energy is independent of momentum (i.e. independent of 
e). The arguments that justify this simplification are 
provided, for example, in Ref. [5|- 

In Eq. 0, if N (e) is taken to be a constant (= 
N o (0)), extending over all energies, then using the fact 
that X is independent of momentum (i.e. e), we obtain 
N(lu) = N o (0). Thus, the standard approximation, that 
the Fermi energy and bandwidth are large energies com- 
pared to the phonon energy (so that we can neglect the 
former and simply integrate from —00 to 00) leads to 
an electron density of states which is unmodified by the 
electron-phonon interaction. This is true even though the 
self energy has a non-trivial frequency dependence fl3| . 
When a more realistic bare electron density of states is 
used, then Eq. leads to an altered EDOS. For exam- 
ple, with N (e) — N o (0)8(D jl — |e|), i.e. a constant over 
a limited energy range, Eq. gives 
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If a Lorentzian form is used the bare density of states is 
given by 



(D/2) 2 



(7) 



3 



In either case -/V o (0) is the density of states at the Fermi 
level and D is the full bandwidth, defined in an obvious 
way in the case of the constant case, and as the full width 
at half maximum in the Lorentzian case. For this latter 
case, Eq. (0 gives 



D/2(D/2 + \TmE(u + id)\ 



3 - 



N(m) _ 

N o {0) ~ (RcS(w + iS) - lu) 2 + (|Im£(w + iS)\ + D/2) 2 



.(8) 
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Eq. © or © is required to self-consistently calculate 
the electron self energy given by Eq. ■ I n what follows 
we further simplify the calculation by adopting T = 0, 
as in Ref. |l2|. This simplifies Eq. since f(ui) — > 
0(—u>) and n(v) — * in this limit. Also, we will adopt 
particle-hole symmetry throughout this paper. The final 
equations, separated out into their real and imaginary 
parts, are given by 
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FIG. 2: The sequence of phonon spectral functions used to 
model broadening. The curves are truncated Lorentzians as 
described in the text with widths as labelled. 
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and the logarithmic singularity is now explicit. Using 
an EDOS that is constant with infinite bandwidth gives 
zero for the lu' integral in Eq. 1|12|) . and we recover the 
'standard' result Q for the electron self energy. 



N(lu + v) 
JV o (0) 



(1 - 0(u + v)) 



As mentioned in the introduction, the primary purpose 
of this paper is to explore the consequences of a non- 
singular phonon spectrum. For simplicity we will adopt 
a Lorentzian spectrum for the phonons, given by 



a 2 F{v) 
\'u>e 



2ir \{v-lu e ) 2 + S 2 



S 2 
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where the subtracted term ensures that the spectrum 
is continuous everywhere, particularly at the endpoints, 
and n is the full width of the spectral function (less than 
loe)- As the parameter 5 approaches zero, this spectrum 
approaches an Einstein spectrum centered at lu = lue, 
and A' — ► A. Fig. 2 shows several spectra for different 
values of 6. 

As a technical aside, simplifications occur for the Ein- 
stein spectrum, where a 2 F{v) = ^^-5(v — u>e)- As is 
apparent from Eq. (|10() . Im£(cj + id) becomes simply 
related to the self-consistent EDOS. The real part of the 
self energy becomes singular, as is evident when Eq. 
is rewritten for an Einstein phonon spectrum as 



Re£(w + id) 



Xujf 



, ,{N(lu')-N{uj-lje)- 
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N(u) — uje) 
JVo(0) 



ln| 



J uj 2 - (uje+uj 1 ) 2 
uj — UJE 



LU + LUE 



(10) RESULTS 

Einstein phonon spectrum 

It is clear from the discussion in the previous section 
that there is a simple scaling relation amongst the ener- 
gies in the problem. Nonetheless wc will use real units, 
and the reader can scale the results to other energy scales, 
if so desired. We begin with lue = 10 meV, and use 
a bandwidth 10X this amount, i.e. D = IQlue = 100 
meV. For dcfinitcness we use A = 2, which is considered 
very strong coupling (Pb, for example, has A f» 1.5), and 
AT o (0) = l/D. In Fig. 3 we plot the (a) real and (b) 
imaginary parts of the electron self energy as a function 
of frequency for these parameters. We adopt a Lorentzian 
shape for the bare EDOS. Three curves are shown; one 
is for the standard theory, where an infinitely wide band 
with constant density of states is assumed, the second is 
for the non-self-consistent result, where the bare EDOS 
(N (lu)) is substituted for N(lu) on the right hand side 
of Eqs. (|12I10|I . and the third represents the full self- 
consistent solution to these same equations, using Eq. 
iJSJl in addition. 

In parts (c) and (d) we plot the same quantities, cal- 
culated this time for a bare EDOS which is constant be- 
tween —D/2 and D/2. What is clear from these plots 
is that a singularity in the real part of the self energy 
exists at the Einstein frequency, regardless of the par- 
ticular approximation used. In fact, as Eq. I|12|) makes 
clear, the singularity is logarithmic, and exists regardless 
of the value of A. Thus, even for A « 0.01, as long as 
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FIG. 3: The (a) real and (b) imaginary parts of the elec- 
tron self energy as a function of frequency, using an Ein- 
stein spectrum for the phonons and a Lorentzian for the bare 
EDOS, with uje/D = 1/10. The curves are computed with 
loe = 10 meV, and A = 2. The non-self-consistent result 
uses the diagram shown in Fig. 1 with the fully interact- 
ing electron Green function replaced with the non-interacting 
electron Green function. Also shown is the 'standard' result 
with infinite electron bandwidth (short-dashed curve). Note 
the singularity in the real part of the self energy, which occurs 
at all levels of approximation. Parts (c) and (d) show similar 
results obtained with a constant bare EDOS with bandwidth 
D. An additional singularity occurs at ui = D + cue and is 
due to the abrupt band-edge in this model. Note that the 
imaginary parts of the self energy, (b) and (d), are zero up 
to uje, and then turn on abruptly, because the phonon is a 8- 
function. For the constant density of states (d) this imaginary 
part of the self energy extends to infinity in the case where the 
bandwidth is infinite (dotted curve). For a finite bandwidth 
(dashed curve) , the imaginary part of the self energy extends 
to D + LOE- 



it is nonzero, the Migdal approximation results in a log- 
arithmic singularity in the self energy. That this is not 
a serious problem is hinted at by Eq. JjJJl, where one 
can see that, as long as a broader phonon spectrum is 
used, the logarithmic singularity will be integrated to a 
non-singular result. Furthermore, as mentioned earlier, 
even in the case of an Einstein phonon spectrum, the 
electron density of states will remain unaltered when the 
'standard' approximation of infinite bandwidth is used. 
The point of Ref. |l2( was, however, that this is not the 
case when a bare EDOS with a non-infinite bandwidth 
is used. The EDOS is plotted in Fig. 4, again using the 
same parameters as in Fig. 3, and for the same levels of 
approximation (the bare EDOS is also included for ref- 
erence). In Fig. 4a (b) we use a bare EDOS which is 
Lorentzian (constant) with bandwidth D. In both cases 
the singularity manifests itself in the final EDOS in both 
the non-self-consistent and self-consistent Migdal approx- 
imations. Thus, it would appear that the EDOS has 
collapsed, and the effective bandwidth is of order 2lve- 
However, an examination of the self-consistently deter- 
mined EDOS is shown in Fig. 5, for several values of 
A. Here the collapse of the band is shown explicitly to 
occur for even very small values of A, consistent with the 
singular behaviour in the real part of the self energy. Yet 
the Migdal approximation ought to be valid at least for 
very weak coupling. 

Fig. 4 shows another interesting feature, which is the 
significant alteration of high energy states. This is partic- 
ularly prominent in Fig. 4b, where states are created at 
energies above the band edge. This would occur for much 
smaller values of the electron-phonon coupling strength, 
and for even larger values of the bare bandwidth, D. 
This result is somewhat counterintuitive. We normally 
anticipate that a perturbative interaction affects states 
just near the Fermi level. However, here, as in the exact 
solution, all states are modified in an additive way, so 
even states well away from the Fermi level get pushed to 
higher energies. 

Rather than focusing on the self energy correction it- 
self, Alexandrov et al. used a different criterion for 
the phenomenon of band collapse (which they attributed 
to polaron formation). They simply took the full width 
at half maximum of the converged EDOS, regardless of 
what structure the EDOS contains at lower frequencies. 
For example, in Fig. 5 there would first be an initial in- 
crease in the effective bandwidth as A increases. Only for 
A > 3 does the effective bandwidth decrease below D (in 

this case 100 meV). For increased coupling strengths, the 
effective bandwidth decreases smoothly to about 4D/5, 
but then has a number of erratic jumps as it decreases 
to approximately 2u>e- Note that Alexandrov et al. fl^ ] 
found a smooth decrease to \loe for two reasons. First, 
much of the fine structure visible in Fig. 5 was not ob- 
tained in their solutions, and second, they included a sec- 
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FIG. 4: The self-consistently computed EDOS for the pa- 
rameters of Fig. 3, with (a) a Lorentzian bare EDOS, and (b) 
a constant bare EDOS. Note the collapse at u> — u)e, and the 
multi-phonon structure apparent in the self-consistent calcu- 
lation. 



ondary band which served to smooth out some of the fine 
structure and reduce the impact of the electron-phonon 
coupling. 

Nonetheless Fig. 5 demonstrates that the Einstein 
model for the phonon spectrum leads to anomalous be- 
haviour in the self-consistent EDOS; to determine how 
much of this is due to the physics of band narrowing, 
and how much is attributable to the singular nature of 
the phonon spectrum, we will study the effect of a broad- 
ened phonon spectrum. 

Before doing so, however, we show in Fig. 6 the self- 
consistently calculated EDOS for several values of Ue/D. 
In Fig. 6a (6b) we use A = 1.0 (5.0), and plot the result- 
ing EDOS for uj e /D = 1/2, 1/5, 1/10, and 1/20. Clearly 
as the Einstein frequency becomes comparable to the 
bare bandwidth the narrowing effects become more pro- 
nounced, particularly for large values of A. Interestingly, 
while in the opposite limit, uje/D — > 0, we approach the 
'standard' model where the bare EDOS is unmodified by 
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FIG. 5: The self-consistently computed EDOS computed 
with a bare Lorentzian EDOS with lue/D = 1/10, com- 
puted for a variety of values of the electron-phonon coupling 
strength. Note that the collapse at u = lue occurs for all val- 
ues of A. The multi-phonon structure becomes more apparent 
as A increases. Moreover, more spectral weight is pushed to 
higher energies as A increases. 



interactions, this particular approach to that limit always 
shows the collapse in the EDOS at u> — loe- This is true 
even in the case of a constant bare band with finite width, 
where there is an even closer connection to the standard 
model in this limit. 
Lorentzian phonon spectrum 

In Fig. 7 we show (a) the real part and (b) the imag- 
inary part of the self energy and (c) the self-consistent 
density of states, for a band of electrons with a bare 
Lorentzian EDOS with width D = \Qloe interacting with 
a broadened electron-phonon spectrum (5 = 5.0 meV 
in the phonon spectrum Lorentzian centred at we = 10 
meV). Results are shown for A = 0,0.1,0.5,1.0,2.0 and 
3.0. By the criterion described above that was used by 
Alexandrov et al. |l2| . no band narrowing has occurred 
up to A = 3.0. However, inspection of the self-consistent 
EDOS shown illustrates that this criterion may be too 
simplistic to describe the more global behaviour that is 
evident in Fig. 7. 

For instance, it is clear that profound changes take 
place within loe oi the Fermi level. So, while the full 
width at half maximum actually increases with increasing 
A in the range shown, the EDOS clearly narrows close the 
Fermi energy (uj — 0). Thus we could plot instead the 
normalized EDOS at lu w loe, for example, as a function 
of coupling strength. This is shown in Fig. 8 for various 
values oi uje/D. 

Interestingly, Fig. 8 shows that the fastest modifica- 
tion with increasing coupling strength occurs near A = 0. 
This is in contrast to the criterion used in Ref . 0] > where 
the measure of band collapse used there plummeted at 
a particular value of A ~ 3. The development of a reso- 
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FIG. 6: The self-consistently computed EDOS computed 
with a bare Lorentzian EDOS with (a) A = 1 and (b) A = 5. 
The latter is chosen purposely very large to better see the 
trends. We use four values of the adiabatic ratio loe/D, as la- 
belled. In (a) the bare Lorentzian EDOS are also shown with 
lighter curves so that the modifications due to the electron- 
phonon interaction are readily visible. The bare EDOS are 
not shown in (b) for clarity; nonetheless, it is clear that more 
significant changes occur for the higher value of A. Note that 
there is not a sum rule when the bandwidth D is changed, as 
is the case here. 



nant peak near the Fermi level can indicate new physics 
(e.g. polaron formation), but there is nothing in Fig. 8 
to indicate that this occurs abruptly at some coupling 
strength. Calculations with a bare electron band which 
is constant with bandwidth D produce results qualita- 
tively and quantitatively similar to those shown in Fig. 
8. This result clearly shows a very smooth evolution of 
the EDOS as a function of coupling strength out to very 
large values of A. Inspection of the self-consistent EDOS 
as a function of frequency for values of A near 10 show 
no qualitative differences with those shown in Fig. 7c. 
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FIG. 7: (a) Real and (b) Imaginary parts of the electron 
self energy vs. frequency for a bare Lorentzian EDOS with 
uje/D = 1/10, using a broadened truncated Lorentzian line- 
shape for the phonon spectrum. We use central frequency 
uje = 10 meV, halfwidth 5 = 5 meV, with truncation at 
r/ = cue ± 9.95 meV. Results are shown for various values of A 
as labelled. Note that a singularity in (a) is no longer present. 
Panel (c) shows the self-consistently computed EDOS using 
the results of (a) and (b). Consistent with (a) a collapse at 
u) = cue is no longer present. All sharp features are no longer 
present, but a dip near to — uje and excess spectral weight at 
high energy remain. 
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FIG. 8: Normalized EDOS at uo = uje vs. coupling strength. 
This frequency is used because it roughly corresponds to the 
minimum in the EDOS (see Fig. 7c). Note that the collapse 
is most severe when the bare bandwidth is comparable to 
the characteristic phonon frequency. Also note that the "col- 
lapse", as measured by this indicator, evolves smoothly as a 
function of coupling strength. The most significant change 
occurs at weak coupling. 



SUMMARY 

We have revisited the Migdal approximation at a 
slightly more sophisticated level than the 'standard' 
treatment where an infinitely broad electron band is as- 
sumed. This was done by self-consistently computing the 
electron density of states for electrons in a band with fi- 
nite bandwidth interacting with phonons. In fact the self- 
consistency is not at all necessary to observe the changes 
near the Fermi level that result from this interaction. The 
non-self-consistent calculation captures the tendency for 
the band to form a resonance with width given by the 
characteristic phonon energy just as well. The key ele- 
ment is that the bare EDOS has a finite bandwidth. 

Previous work nas focused on the Einstein spec- 
trum for the phonons coupled to an electron band de- 
scribed by a Lorentzian density of states. We have con- 
sidered a constant bare EDOS as well, and found very 
little difference in the results. A qualitative change in the 
results does occur, however, when a broad phonon spec- 
trum is used instead of the delta function that character- 
izes the Einstein model. In the latter case the electron self 
energy is always singular, regardless of the level of self- 
consistency used to calculate the Migdal approximation. 
With a bare EDOS with finite bandwidth, this singular- 
ity results in an electron density of states that collapses 
at the Einstein frequency loe, for any nonzero value of 
the electron-phonon coupling constant, A. This collapse, 
however, is due to the unphysical nature of the Einstein 
spectrum [Tij , and does not signal a metal-insulator tran- 
sition. 



As expected, the use of a broadened phonon spectrum 
eliminates the singularity in the self energy, and in the 
self-consistent EDOS. Fig. 7c epitomizes the change that 
occurs (compare with Fig. 5). There is still a strong sup- 
pression of the EDOS at energies u>e away from the Fermi 
energy, which indicates that a resonance occurs near the 
Fermi level. These states clearly are pushed to much 
higher energies. The fact that the energy scale for the 
resonance is u>e is indicative of increased involvement of 
phonons in the electron states near the Fermi level (and 
therefore of polaron formation) , but there is no signal in 
these results of a collapse of the conduction band for A 
of the order of unity. Eventually, as Fig. 8 indicates, the 
resonance dominates the EDOS for very large values of 
A, particularly for relatively large values of the adiabatic 
ratio, oje/D. In fact, for small values of the adiabatic 
ratio, the self-consistent Migdal approximation is prop- 
erly adjusting the EDOS to at least partially incorporate 
some of the physics of polaron formation, i.e. that the en- 
ergy scale for the electrons becomes comparable to that 
of the phonons. 

Ref. [ijj utilized a secondary band; as the presence 
of this band serves to 'soften' the impact of the electron- 
phonon interaction on the primary band, we have omitted 
it here. The same qualitative results are obtained, except 
for somewhat higher values of the coupling strength [15( . 

Finally, over the last four decades there have been 
many studies of electrons interacting with phonons, and 
the potential of a crossover to a regime where polaronic 
behaviour dominates the physics. Many of these stud- 
ies use the Holstein model for the phonons, for the sake 
of simplicity. The exact studies (see Ref. for a 

short review and pertinent references) account for the 
phonon renormalization; thus, in principle, these include 
phonon broadening effects. In practice, unfortunately, 
many of these studies are carried out on finite lattices, 
or in various limits (e.g. the adiabatic approximation), 
so that a completely satisfactory solution is not avail- 
able. Nonetheless, as reviewed in Ref. the majority 
of these studies suggest that a cross-over occurs from 
free electron-like to polaronic behaviour, near A ~ 1. 
We make a cautionary note, however, in relation to this 
work; the previous statement applies to the bare dimen- 
sionless coupling constant, A. However, in those studies, 
the 'operational' value of A is in fact much higher, be- 
cause phonons have softened, etc. 0. It is this 'opera- 
tional' value which more closely corresponds to the value 
of A used in this work. In any event, the limitations 
on coupling strength in the Migdal approximation in the 
normal state, or the corresponding Eliashberg formalism 
in the superconducting state [2(| must ultimately come 
from these exact studies. Achieving a non-singular result 
within the Migdal approximation scheme is insufficient 
grounds for its accuracy. 

We should also remark that some work has also been 
done on the Barisic-Labbe-Friedel model ^} (also known 
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as the SSH (Su-Schrieffer-Heeger) model H3)j where dis- 
persion in the phonon spectrum exists at the start. This 
model may be distinct from other cases where phonons 
are broadened due to the electron-phonon interaction 
itself or anharmonicity, since it contains sharp disper- 
sive modes vs. individually broadened dispersive modes. 
Questions concerning these models deserve further study. 

This work was supported by the Natural Sciences and 
Engineering Research Council (NSERC) of Canada and 
the Canadian Institute for Advanced Research (CIAR). 



[1] A.B. Migdal, Zh. Eksp. Teor. Fiz. 34, 1438 (1958) [Sov. 

Phys. JETP 7, 996 (1958)]. 
[2] After all, the vertex correction should be compared with 

a diagram that is of the same order in the coupling and 

is retained. 

[3] S. Engelsberg and J.R. Schrieffer, Phys. Rev. 131 993 
(1963). 

[4] We use the term 'Migdal approximation', although in the 
literature this is often referred to as 'Migdal's Theorem'. 

[5] For a review, see, for example, P.B. Allen and B. 
Mitrovic, in Solid State Physics, edited by H. Ehrenreich, 
F. Seitz, and D. Turnbull (Academic, New York, 1982) 
Vol. 37, p.l. 

[6] For a text, see, for example, S. Nakajima, Y. Toy- 
ozawa, and R. Abe, The Physics of Elementary Exci- 
tations (Springer- Verlag, New York, 1980) . 

[7] M. Abramowitz and LA. Stegun, Handbook of Mathemat- 
ical Functions (Dover, New York, 1964). 

[8] D.J. Scalapino, in Superconductivity, edited by R.D. 
Parks (Marcel Dekker, Inc., New York, 1969)p. 449. 

[9] B.N. Brockhouse, T. Arase, G. Caglioti, K.R. Rao and 

A. D.B. Woods, Phys. Rev. 128 1099 (1962). 

[10] Migdal Q began with a bare coupling strength Ao , which 
led to strong renormalizations of both the phonon and 
electron properties. In following Ref. fl^l we compute 
electron properties only, so that our coupling strength 
is already renormalized in some sense, to the value that 
Migdal calls A. While Ao was found to be required to be 
less than unity to preserve phonon properties, A, can, in 
principle, take on larger values. 

[11] For a nice exposition of the formalism, see, for example, 

B. Mitrovic and J. P. Carbotte, Can. J. Phys. 61, 758 
(1983). 

[12] A.S. Alexandrov, V.N. Grebenev, and E.A. Mazur, 
Pis'ma Zh. Eksp. Teor. Fiz. 45 357 (1987) [JETP Lett. 
45 455 (1987)]. 

[13] This statement applies to a calculation of the electron 
density of states. However, it is worth mentioning that a 
calculation of the electronic specific heat at low tempera- 
tures, for example, leads to the usual expression for a free 
electron gas, modified by the mass enhancement factor, 
1 + A. Thus, as far as the specific heat is concerned it looks 
as if the electron density of states has been modified by 
the electron phonon interaction. This occurs in this and 
other properties because the calculations involve more 
than just replacing the free electron density of states by 
the one modified by the electron phonon interaction. 

[14] A similar 'collapse' occurs in the free electron density 



of states, computed in the Hartree-Fock approximation. 
This is remedied by invoking screening. See, for exam- 
ple, N.W. Ashcroft and N.D. Mermin Solid State Physics 
(Saunders College Publishing, New York 1976), particu- 
larly the discussion on p. 344. 

[15] F. Dogan, Masters Thesis, University of Alberta, Depart- 
ment of Physics (2002). 

[16] A.S. Alexandrov, Europhys. Lett. 56, 92 (2001). 

[17] S. Barisic, J. Labbe and J. Friedel, Phys. Rev. Lett. 25, 
919 (1970). 

[18] W.-P. Su, J.R. Schrieffer, and A.J. Heeger, Phys. Rev. 
Lett. 42 1698 (1979); Phys. Rev. B22, 2099 (1980). 

[19] F. Marsiglio, Phys. Rev. B 42 2416 (1990). 

[20] For example, some authors have investigated the so- 
called asymptotic limit of Eliashberg theory, for large 
values of A. See, for example, P.B. Allen and R.C. Dynes, 
Phys. Rev. B12 905 (1975), JR. Carbotte, F. Marsiglio 
and B. Mitrovic, Phys. Rev. B33 6135 (1986), and R. 
Combescot, O. V. Dolgov, D. Rainer, and S. V. Shulga, 
Phys. Rev. B53 2739 (1996). These and many other re- 
sults provide no justification for the use of Eliashberg 
theory. Justification, or demonstration of the invalidity 
of Eliashberg theory, as in the case of the Migdal ap- 
proximation, must come from more exact studies. 



